Metastable states and macroscopic quantum tunneling in a cold atom Josephson ring 
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We study macroscopic properties of a system of weakly interacting neutral bosons confined in 
a ring-shaped potential with a Josephson junction. We derive an effective low energy action for 
this system and evaluate its properties. In particular we find that the system possesses a set of 
metastable current-carrying states and evaluate the rates of transitions between these states due to 
macroscopic quantum tunneling. Finally we discuss signatures of different metastable states in the 
time-of-flight images and argue that the effect is observable within currently available experimental 
technique. 
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Macroscopic Quantum Tunneling (MQT) is undoubt- 
edly one of the most remarkable manifestations of macro- 
scopic quantum coherence, a property intrinsic to su- 
perconductors and Bose-Einstein Condensate (BEC) sys- 
tems. The most well studied system exhibiting MQT is 
the so-called phase qubit, a superconducting loop con- 
taining a Josephson junction, where MQT occurs be- 
tween metastable current-carrying states of the Cooper 
pair condensate in the loop [l|, [2[. The main mecha- 
nisms controlling these metastable states and the MQT 
arise from an interplay between several factors, such as 
Josephson, inductive, and charging energies, as well as 
dissipation due to shunting resistances and/or quasipar- 
ticle excitations. The Josephson and the inductive ener- 
gies form an effective potential for the flux-phase variable 
(i.e., the phase difference across the Josephson junction), 
while the charging energy (resulting from plasma oscilla- 
tions due to the Coulomb interactions between electrons) 
gives rise to the effective kinetic energy for the phase 
[2, Q. Dissipation is also known to play a critical role 
in MQT dynamics, leading to a substantial slow-down in 
tunneling rates [4|, [5[ as well as to the localization of the 
magnetic flux in the case of symmetric double well effec- 
tive potential (realizable by applying a 1/2 magnetic flux 
quantum through the loop) [g, Lfl] • 

In this Letter we consider dynamics and properties of 
a phase qubit-like (Josephson ring) device based on a 
BEC of neutral particles. Contrary to the conventional 
superconductor-based devices, these systems did not re- 
ceive much attention in literature, - presumably due to 
lack of their experimental realizations in the past. Our 
interest in such systems, however, is not purely academic: 
advances in cold atom trapping, such as the recently 
demonstrated "painted potential" technique [§], open a 
clear path for the realization and development of such 
phase qubits. As we will show below, these cold atom per- 
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FIG. 1: (Color online) (a) The system of interest: Density 
profile of a BEC in a "painted potential", i.e., a ring with 
a tunnel junction (see Ref. 0] for comparison with actual 
experimental photographs), (b) and (c) Effective potentials 
as a functions of phases across the Josephson junctions for 
two cases: (b) Q = 0, a = 5.8 and (c) Q ^ 0, a = 5.8 and 
6o = 5.8. 



sistent current devices posses a number of properties sim- 
ilar to their superconducting counterparts. These prop- 
erties, however, are controlled by rather different physi- 
cal mechanisms, which rely not on particles' charge (like 
inductive and charging energies of the in superconduct- 
ing devices), but on the effects of quantization due to 
the finite size of the system. Indeed, while for super- 
conductors the phase (Goldstone) mode is gapped due to 
plasmon oscillations, it is gapless for the neutral particles 
interacting through a short range potential. As a result, 
the phase dynamics for the latter systems is susceptible 
to the long-range fluctuations and is strongly affected by 
the system's dimensionality, geometry of the traps, etc. 

We consider a BEC confined in a ring-shaped poten- 
tial containing a single cut, i.e., a tunnel (Josephson) 
junction, see Fig. 1(a). We derive an effective action 
for the phase difference across the Josephson junction. 
Dynamics described by this action strongly depends on 
the system's size (i.e., the ring's circumference), and is 
dramatically different in the two limiting cases, i.e., for 
small and large BEC rings. In the latter case dynamics 
of the phase <p across the junction corresponds to that of 
an overdamped particle, diffusing in a periodic potential 
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(due to the Josephson energy at the junction, see below). 
In the former case the "phase-particle" is underdamped 
and moves in the potential, which is a sum of Josephson 
energy (~ cos(/>) and effective inductive energy (~ <j) 2 ). 
The inductive energy term is of quantum origin: It arises 
due the kinetic energy of the BEC particles confined in 
the ring. The competition between Josephson and the 
effective inductive energies gives rise to the formation 
of metastable current-carrying states, i.e., with 0^0. 
We obtain criteria for the formation of these states and 
discuss means of their control and possibility of their re- 
alization in cold atom systems. We analyze MQT rates 
from these states as a function of system's parameters for 
stationary and for rotating potentials. Finally we out- 
line signatures of these current-carrying states in time- 
of- flight experiments. 

The Hamiltonian for the system of interest is given by 
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where potential V(x) describes the tunnel barrier and the 
system is assumed to be one dimensional [9] with periodic 
boundary conditions. Here and in the following we use 
units with h = 1 and fc# = 1, unless stated otherwise. 
We would like to derive an effective low energy action 
for the system described by the Hamiltonian in Eq. (pQ). 
This can be carried out as follows. In the absence of 
the potential barrier V(x) this action is readily obtained 
by representing Bose field operators in terms of hydro- 
dynamic variables as ip = \J p + 5pe l(p , where p = N/L 
and N is total number of atoms in the ring. In that case 
the linearization with respect to the small density vari- 
ations 5p yields an effective phonon Lagrangian in the 
long wavelength limit [lQj |. In the phase representation 
it reads 
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where c is sound velocity for an equivalent infinite sys- 
tem (c = (pA/m) 1 / 2 in the weak coupling limit). Here 
and in the following we use imaginary (Matsubara) time 
representation. The spectrum of the phonons is discrete 
due to the finite size of the system, with energy spacing 
between the adjacent levels AE = 2irc/L. If the barrier 
is sufficiently high, the inclusion of the tunnel barrier po- 
tential V(x) changes the functional dependence of the La- 
grangian density on ip only locally. That is, the BEC La- 
grangian density in the region with V{x) =0 is still given 
by Eq. ([2j), supplemented by the appropriate boundary 
conditions at the barrier (see below). For the high poten- 
tial barrier transitions across the barrier can be described 
by the usual hopping term ~ [ip^ (O)ijj(L) + ^ (L)ip(0)] 
and therefore the effective action can be cast in the form 
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where <j) = cp(0) — ip(L) is the phase difference across 
the junction and t is the transmission coefficient. Obvi- 
ously the value of t depends on the shape of the poten- 
tial barrier as well as it accounts for the depletion of the 
density near the barrier. Moreover, it may depend on 
the interactions between particles and their density [Tlj . 
However, if the height of the barrier Vb is large com- 
pared to the interaction energy Ap, this dependence can 
be neglected, and, as a rough estimate, one can use the 
standard semiclassical expression for the rectangular bar- 
rier, t ~ (1/md) exp [— d(2mVo) 1//2 ], where d is the width 
of the barrier. Furthermore, Eqs. ([21 [3|) are applicable 
only when Josephson energy per particle pt/N is small 
compared to the level spacing AE, i.e., for c^> t; other- 
wise the presence of the barrier significantly affects the 
quasiparticle (i.e. phonon) states inside the ring and the 
weak link picture is no longer valid. In this weak link 
limit the supercurrent at x = and x = L is vanishingly 
small, and therefore Eq. ([3]) must be supplemented by 
the boundary conditions d x cp(0 : r) = d x cp(L,r) = 0. 

The action given by Eqs. ([21 [3]) can be reduced to a 
local action by integrating out the <p(x,r) field every- 
where except x — and x = L. This is accomplished by 
introducing the functional (5-function according to 

Z = n / ^^e- 5 (5[^)-(^(0,^) + (^(L,cj)], (4) 

and using identity 5(a) = (2tt)~ 1 J dz exp (iza). The 
fields cp and z can then be integrated out (note that in 
doing so one should expand field ip(x,r) in cos(7rnx/L), 
thus satisfying the boundary conditions above). After 
some calculation we obtain that the partition function 
takes the form Z = J P^exp (— iS e ff), where 
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dr cos (j)(r). (5) 



Effective action in Eq. ([5]) is similar to that for an 
impurity in the Luttinger liquid [12[. The first term of 
the right-hand side is different due to the finite length of 
the ring. Indeed, for sufficiently large ring the first term is 
proportional to \u\ and the action is that of a dissipative 
particle in a periodic potential. In the opposite small 
ring limit, the first term in Eq. ([5]) yields 
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where the effective kinetic energy term describes the low- 
est BEC mode, that is, motion of the BEC as a whole, 
and the last "inductive" term arises due to the kinetic 
energy of the BEC particles in the ring. Together with 
the cosine term in Eq. ([5]) this <j) 2 term provides effec- 
tive potential energy for the "phase particle". This po- 
tential may have one (at (j) = 0) or more (at <p ^ 0) 
local minima. The <p ^ minima correspond to the 
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metastable states carrying non-zero current. The shape 
of such effective potential can be conveniently character- 
ized in terms of the dimensionless parameter a = mtL. 
For a — ot\ ~ 4.60 the first two metastable minima ap- 
pear at (j) « ±1.437r. Two more metastable states appear 
when a > <%2 ~ 10.95, and so on. Furthermore, the effec- 
tive potential can be modified by rotating the "painted 
potential" . As a result the potential energy of the super- 
fluid in the rotating reference frame changes by — MO, 
where O is angular frequency of rotation and M is angu- 
lar momentum of the superfluid (in the stationary frame) 
Pij, and therefore the effective potential V(4>) can be cast 
in the form 
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where </>o = 2irIQ (here I is the moment of inertia of 
the superfluid per atom, / ~ mL 2 ). Such V eE ((j)) is 
equivalent to that of a Josephson flux qubit, i.e., a su- 
perconducting loop with a Josephson junction in a mag- 
netic field [7]. For </>o = 7r effective potential V eff (</>) 
is a symmetric double well, while for <j)$ > tt the mini- 
mum close to becomes metastable and then, upon fur- 
ther increase of O, disappears completely for 0o > </>o r = 
[ Va 2 -1 + 1 cos _1 ( -l/a) \}/a. 

Note that in Eq. (JJ) the particle density controls only 
the overall strength of the effective potential, but not 
its relative shape. The latter depends only on the pa- 
rameter a, which is determined by the properties of the 
"painted potential", but not by the number of particles 
in the system. This fact is rather advantageous from 
the experimental point of view: While the number of 
particles in the system typically obeys Poissonian statis- 
tics and therefore it can significantly fluctuate in differ- 
ent experimental runs, the form of the painted potential 
(i.e., the value of parameter a determining the onset of 
the metastability) is easy to control with sufficiently high 
precision. 

Transitions between different minima of the effective 
potential in Eq. Q are driven by thermal and quantum 
fluctuations. The former mechanism dominates at suffi- 
ciently high temperatures, exceeding frequency of small 
Josephson oscillations in a given well. Note that in the 
metastable region, i.e. for mtL > 1, this frequency is 
of the order of the spacing between quasiparticle energy 
levels, AE. For the system of 87 Rb atoms L ^ 0.2 mm 
and p 3D ~ 10 12 cm- 3 , e.g. Ref. [8], T* ~ AE ~ 0.2 nK. 
Thus, above T* the rate of transitions is determined by 
the thermal activation rate ~ exp (— F eff /T), where Vq S 
is the height of the effective potential barrier in Eq. ([7j). 
In the following we evaluate transition rates for two cases: 
(I) For a system in the first metastable state, i.e., for 
ol\ < a < a 2 and O = 0, e.g. Fig. 1(b); (II) For the 
rotating system (Q ^ 0), when the initially global mini- 
mum at (j) = becomes metastable, i.e., for tt < <po < ^q 1 *, 
e.g., Fig. 1(c). In the former case the rate of thermally 
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while for the latter case 
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For the order of magnitude estimate both pre-exponential 
factors r and T f can be taken as frequency of small 
Josephson oscillations ~ AE [l3j], while the height of 
the effective potential barrier V r eff ~ p/(mL) in Eqs. ([H 
[9]) is of the order of 10 nK (here and in the following 
we use the same values for the system's parameters as 
above; see also Ref. [8[). Thus we expect that thermally 
activated transition rates of the order of 1 — 10 Hz are 
readily observable at temperatures below 10 nK. 

For temperatures below T* ~ 0.2 nK thermal ac- 
tivation mechanism becomes ineffective and transitions 
are driven by quantum fluctuations, i.e., by the MQT. 
The rate of the MQT can be computed by evaluating 
Euclidean action over the classical instanton trajectory 
T = r o exp[-5(0 c /)] [14]. Measurable MQT rates can 
be found near the spinodal instabilities of the metastable 
states, i.e., when the effective potential barriers are rela- 
tively small. In this regime one can expand the effective 
potentials around the local minima (state 2 in Fig. 1(b), 
e.g. case (I), and state 3 in Fig. 1(c), case (II)) up to 
the cubic terms in (</> — </> m in)- The kinetic energy, i.e., 
the first term in Eq. ([5]) minus its value at uo = (the 
second term in Eq. is non-local in Matsubara time 
r and therefore is difficult to be dealt with exactly. For, 
ql\ — a <C 1, however, the instanton trajectory (j>ci(j) 
contains only small frequencies. In this limit the kinetic 
energy is local, e.g., the first term in Eq. (|6]). In the 
non-local ql\ — a > 1 regime we use variational approach 
by choosing (j) c i{j) = (/>% ex P ( — I t I/At), where (jP cl and 
At are variational parameters. The results for O = 0, 
i.e., case (I), are 

In El = n J 2 -° 7(a " ai ^ /A a ~ ai <<C 1 h m 
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where r] = y pj (mX) and the pre-exponential factor To 
is again of the order of AE [15[ . Evaluation of the MQT 
for the rotating system, i.e., case (II), gives 
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Note that the exponents for the MQT rates in Eqs. ([TOl 
[TT]) are not proportional to the system's size or its total 
number of particles, e.g. Ref. [16[. Indeed, while the 
effective mass, e.g. Eq. (|6j), is proportional to L, the 
height of the effective potential barrier scales as 
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In summary, we have proposed a novel BEC based 
qubit-like device and theoretically studied its properties. 
We have identified a set of macroscopic metastable states 
that such system is expected to exhibit under certain 
conditions and evaluated transition rates between these 
states. We have discussed signatures of these states in 
the standard TOF measurements and argue that the ef- 
fects considered in this paper can be readily observed in 
contemporary cold atom systems. 

We thank I. Martin, V. Privman and E. Timmermans 
for valuable discussions and comments. The work is sup- 
ported by the US DOE. 



FIG. 2: (Color online) Calculated TOF measurement out- 
comes for: (i) Transitions from the metastable current state 
2 — ► 1, Fig. [T] with Q = 0. Observation of (a) corresponds 
to state "1" and (b) — to state "2". (ii) Transitions to the 
current-carrying state in (Q / 0) case. Image (c) corresponds 
to state "3" and (d) — to state "4" . 

The value of parameter r] in Eqs. ([TOf [TT]) for the ex- 
periments in Ref. [8] is ~ 300 and therefore fine tuning of 
parameters is needed in order to obtain reasonable MQT 
rates. While this fact presents a formidable difficulty in 
some MQT proposals, e.g. Ref. [161 ] . it seems that such 
fine tuning is relatively simple for the present system. In- 
deed, the value of parameter 0o is determined solely by 
the system's geometry and by the frequency of rotation, 
but not by the density, and therefore can be controlled 
with desired precision. Moreover, the value of r] can be 
reduced by decreasing the cross-section of the ring [l7| , 
or by increasing the scattering length, i.e., the value of 
A, via the Feshbach resonance. 

Finally we briefly discuss the signatures of the current- 
carrying metastable states in the Time-Of-Flight (TOF) 
measurements. As well known, after the trapping po- 
tential is turned off, at sufficiently long times of expan- 
sion the resulting BEC state corresponds to the Fourier 
transform of the initial state in the trap. In Fig. 2 we 
present numerical Fourier transforms of the initial states 
that have same particle density \i/Jo(r, 0)\ 2 shown in Fig. 
1(a), but multiplied by different phase factors exp (m#), 
where phases 2i\n correspond to the values of (j) in minima 
1 through 4 of the effective potentials in Fig. 1(b) and 
Fig. 1(c). That is, n = for the state 1, while n ~ 0.75 
for the state 2, and so on. The main distinction between 
Fourier images of states with n ~ and with n ~ 1 is 
clear: While the former case corresponds to the momen- 
tum distributions mainly concentrated at the origin, e.g. 
Figs. 2(a) and 2(c), in the later case the distribution is 
ring-like, e.g. Figs. 2(b) and 2(d) [18j . 
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